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NOMENCLATURE 
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D 
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X x »^z 


xz 

J 

K v 

£,JV 

L,N 

N H 

P 

q 


aspect ratio 
lateral acceleration 
state derivative vector 
wing span 

jet-momentum coefficient 
lift coefficient 

lift coefficient of airplane wing 
rolling-moment coefficient 
yawing-moment coefficient 
thrust coefficient 
side-force coefficient 
side-force lag term, 1 - 

factor p sb . g 

lactor (4W) 
acceleration of gravity 

moments of inertia about x- and z-axis 

product of inertia 

cost function (eq. (D24)) 

cold thrust 

side slip vane correction factor 

rolling- and yawing-moment parameters 

rolling- and yawing-moments, transformed variables 

engine speed , % maximum rpm 

roll rate 

dynamic pressure 

ratio of jet-flap chord to wing chord 
yaw rate 
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s 


wing area 


SAS stability augmentation system 

t time 

u,v,w x- 5 y-, and z-components of velocity 
u(t) input vector 

V true airspeed 

v(t) measurement noise vector 
W gross weight 

w(t) process noise vector 

x(t) state vector 

side force 

Y side force, transformed variable 

y(t) output vector 

a angle of attack 

3 angle of sideslip 

6 a aileron deflection 

choke deflection 
6^ landing-flap deflection 

6^ rudder deflection 

6 g spoiler deflection 

6 W control-wheel deflection 

6 wa control-wheel deflection for aileron input 

6^ control-wheel deflection for spoiler and choke input 

6 time-derivative of 6 W ^ 

8 pitch angle 

p ambient air density 

t time; variable of integration 
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<j> roll angle 

Subscripts : 

J due to jet-flap blowing 

o without jet-flap blowing 

p roll rate 

R rudder deflection 

r yaw rate 

w wing, also control-wheel angle 

3 sideslip angle 

5 state and control variable 

Superscript : 

T transpose of vector 

Aerodynamic derivatives (referenced to a system of body axes with the origin at the 
aircraft center of gravity). 
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THE APPLICATION OF PARAMETER ESTIMATION TO FLIGHT MEASUREMENTS 


TO OBTAIN LATERAL-DIRECTIONAL STABILITY DERIVATIVES 
OF AN AUGMENTED JET-FLAP STOL AIRPLANE 
Jack D. Stephenson 
Ames Research Center 


SUMMARY 


Tests were conducted to obtain data for calculating the lateral and directional 
stability derivatives of the Augmented Jet Flap Research Aircraft, an experimental 
short takeoff and landing (STOL) transport vehicle. The derivatives were calculated 
using a linear-regression parameter-estimation procedure, with equations of motion of 
the aircraft represented by a linear model. For the tests reported here, the jet 
flaps of the airplane were deflected 65°, a typical landing approach configuration. 
The study included determination of the effects on the derivatives of angle-of-attack 
variation, the effect of variation in the momentum coefficient of the jet-flap flow, 
and the effect of varying the angular position of the primary jet-thrust nozzles 
(located on the engine nacelles) . This latter variable could be investigated because 
the aircraft is designed so that the pilot can control the nozzle angles to vary the 
thrust direction. 


The stability derivatives are compared with predicted values that had been incor- 
porated into a simulation model for the airplane. At all test conditions, the roll 
damping from the flight measurements was significantly greater than that predicted. 
Effects on the stability derivatives of varying the momentum coefficient of the jet- 
flap flow and of varying the primary jet-thrust angle were relatively small. Values 
obtained for the derivatives and C n ^ were generally in agreement with the pre- 

dictions. The derivatives C n ^ and C n ^ agreed with predictions at moderate angles of 


attack, but their variation with angle of attack differs somewhat from the predic- 
tions. Application of a maximum- likelihood parameter-estimation procedure to the 
flight data yielded average derivative values that are essentially the same as those 
from the linear-regression solutions; however, for some test conditions there was 
considerable dispersion in the maximum- likelihood estimates. 


INTRODUCTION 


The Augmented Jet Flap Research Aircraft was designed and built to provide infor- 
mation on various flight characteristics of a transport aircraft capable of short 
takeoff and landing (STOL) operations. The flight-test program included maneuvers 
that would allow the calculation of lateral-directional stability derivatives for the 
airplane under conditions representing low-speed landing approaches, where, for this 
type of airplane, handling characteristics are generally inferior to those of conven- 
tional aircraft. 

Before the airplane was first flown, its flight characteristics had been pre- 
dicted for use in a simulation model for the Ames Flight Simulator for Advanced 



Aircraft (FSAA) - This model was based on data from tests of one design version in 
the Ames 40- by 80-Foot Wind Tunnel and on predictions from theory of some of the 
stability derivatives. One objective of the present study was to determine the 
accuracy of such predictions as applied to a configuration employing the augmented 
jet-flap, powered-lift concept. The study also demonstrates the use of a linear- 
regression parameter-estimation procedure with a linearized mathematical model for 
the lateral-directional dynamic behavior of the aircraft; the procedure was used to 
obtain estimates of stability and control derivatives throughout a range of angles of 
attack and jet-momentum coefficients. Results obtained by processing some of the 
data using a maximum- likelihood parameter-estimation procedure are presented, and 
these results, as they compare with those from the regression solutions are discussed. 


THE RESEARCH AIRCRAFT 


The Augmented Jet Flap Research Aircraft is a modified high-wing transport 
powered with two nacelle-mounted jet engines. The engines, in addition to providing 
direct propulsive thrust, supply bypass air for a system of lift-augmentation jet 
flaps. This augmentation system is described in references 1 and 2. The engine hot 
exhaust efflux (providing jet thrust) is directed through nozzles (on the sides of 
the nacelles) that can be rotated through an angular range from a few degrees below 
the longitudinal axis of the airplane downward to a position 104° from this axis. 

A photograph and a three-view sketch of the airplane are shown in figures 1 
and 2, respectively. Figure 3 illustrates some of the details of the flap system and 
of the lateral control surfaces. The leading-edge slats (fig. 2) are fixed in the 
extended position. Additional geometric information, weights, and moments of inertia 
are listed in table 1. A more complete description of the aircraft is presented in 
references 1 and 3. 

The lateral-control system utilizes three aerodynamic surfaces: ailerons with 

boundary-layer blowing; spoilers located ahead of the ailerons; and in the outer 
section of each augmentor flap a hinged lower flap element (referred to in this report 
as the "choke"). When deflected, the choke decreases the lift due to the flap. The 
lateral-control surfaces are operated through the control wheel by a central dual 
hydraulic power actuator; the actuator drives one cable system to the ailerons and 
another to the spoilers and chokes. 

Directional control is provided by a two-panel rudder, the aft panel being 
hinged to the trailing edge of the forward panel and geared to it with a 2:1 deflec- 
tion ratio. The rudder is fully powered through a hydraulic actuator controlled by 
cables from the pedals. 

In normal operation the aircraft makes use of a lateral-directional stability 
augmentation system, but for the tests that are the subject of this report, this 
system was switched off. 

The airplane was fitted with a nose boom on which were mounted an angle-of- 
attack vane, a sideslip vane, and Pitot-static instrumentation. An on-board pulse- 
code modulation (PCM) data acquisition system provided time-histories of approxi- 
mately 95 variables at a sample rate of 100 frames per second. In addition to the 
nose-boom measurements, the recorded data included the following: roll, yaw, and 
pitch rates; roll and pitch attitude; linear accelerations in three body axes; 
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ambient conditions; the position of each control surface; engine-performance param- 
eters; and temperatures and pressures in the ducts carrying bypass air to the flap 
system. 




REDUCTION OF DATA 


Aerodynamic interference effects were estimated and corrections were applied to 
airspeed, angle-of-attack, and sideslip data. The airspeed and static pressure 
errors were determined from tests, using a "pacer" helicopter towing a trailing 
"airspeed bomb"; the errors are functions of landing-flap angle and aircraft lift 
coefficient. The angle-of-attack vane error was determined from measurements in 
steady flight of pitch attitude, longitudinal acceleration, airspeed and altitude 
(ref. 1). A correction for the induced sidewash due to sideslip was applied to the 
measurement of sideslip by the nose-boom vane. It was calculated using the equations 
given in appendix A and data from the maneuvers that are the subject of this report. 
The actual sideslip angle was determined to be about 18% smaller than the indicated 
value at angles of attack up to 11°. At higher angles of attack, the correction was 
slightly larger. 

Lateral-control derivatives presented here are derivatives with respect to the 
control-wheel angle, a variable that represents the combined effects of the three 
lateral-control surfaces. The wheel angle (in radians) is used in defining these 
derivatives to permit a direct comparison of the data with predicted control effec- 
tiveness used in the simulation model. However, because of the effects of cable 
stretch and of lags in the actuation system, measured wheel angle was found to be 
unsuitable as a variable to correlate the aerodynamic characteristics of the lateral 
control. For this reason, an effective control-wheel angle was defined based on the 
measured displacements of the control surfaces assuming a linear variation of the 
wheel angle with displacement. Separate wheel angles are computed for each control. 
For the ailerons, defining 6 a as the average of the left and right aileron 
deflection , 


J W; 


3.62 6 a 


Spoilers: 


^ W SP " 1,0 ’"’spoiler 


Chokes : 


6 w C H 3>1 6 choke 


where ^spoiler i s the angular deflection (in degrees) of the spoiler that is 
deployed; is the percent closure of the area between the augmentor flap ele- 

ments; and deflections producing a roll to the left are defined as negative in sign. 
Figure 4 shows the variations of the lateral-control surface deflections with control- 
wheel angle (in degrees), measured with no aerodynamic load and for comparison, the 
effective wheel angles given by the above relations. 

The directional control deflections are computed using the rudder angle, mea- 
sured in degrees, on the forward rudder panel. 

The efflux of engine bypass air from nozzles at the augmented jet flaps produces 
aerodynamic effects that are correlated on the basis of a jet-momentum coefficient 
Cj. This coefficient — equal to J c /(qS), where J c is the ("cold") thrust due to 
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the ejected air — is computed from engine static test data and flight measurements of 
engine rpm, ambient air conditions, and internal duct temperatures and pressures. A 
jet-momentum coefficient, based on the thrust calculated assuming isentropic expan- 
sion of the bypass air flow to static ambient pressure, was determined to be about 
32% larger than Cj defined above at all conditions for which results are presented. 

The engine direct ("hot") thrust (T) and a thrust coefficient C ^ (C^ = T/(qS)) 
were also obtained from engine static test data and from flight-measured engine and 
environmental information. Figure 5 shows, as functions of engine percent rpm, the 
thrust coefficient and the ratio Cj/C^ for an equivalent airspeed of 70 knots and 
standard atmospheric conditions, at two altitudes, 610 M (2,000 ft) and 1,830 m 
(6,000 ft). 


FLIGHT EXPERIMENTS 


Data analyzed in this report were obtained from maneuvers produced by control- 
wheel and rudder inputs approximating doublet functions. During the maneuvers, the 
pilot attempted to maintain either a constant airspeed or a constant angle of attack. 

To perform the maneuver, he established the specified initial conditions in wings- 
level steady flight, applied the control input, and then, with the rudder and lateral 
controls centered, allowed the motion to continue through two to three cycles of 
Dutch roll. 

The flight experiments were conducted at various times as scheduling permitted 
throughout an extended period of time when the aircraft was also being flown for a 
variety of other research purposes. Results of the earlier tests revealed some details 
of the execution of the maneuvers that affected the quality of the data. For example, 
with maneuvers produced by a control-wheel input, a large maximum angular displacement 
of the wheel (greater than about 25°) yielded better quality data than maneuvers that 
were performed with small wheel excursions. In the same type of maneuver, even a 
small movement of the rudder had significant adverse effects. 

All of the tests were performed with the lateral-directional stability augmenta- 
tion system switched off; and without stability augmentation, the handling character- 
istics of the aircraft are poor, particularly at the lowest airspeeds used in the 
tests. As a result, performance of the tests conducted later in the program improved 
as the pilots became more familiar with the flight characteristics of the 
aircraft. In addition to maintaining approximately constant speed, a satisfactory 
maneuver usually required that the Dutch roll response be relatively symmetrical 
about a wings-level attitude, avoiding a tendency toward roll divergence. 

In all of the tests from which data were obtained in this program, the augmented 
jet flaps were set at a deflection of 65°, a typical setting for landing approach. 

With this flap setting the ailerons have 30° droop. Tests were planned to cover a 
large angle-of-attack range below stall and a range of engine thrust levels up to the 
maximum engine rpm specified for continuous operation. To determine the effects of 
thrust vectoring and interactions of the jet efflux and the external airstream, tests 
were conducted with the thrust nozzles at three settings: their position for takeoff 

and cruise, that is 8° downward from the fuselage axis; fully down at an angle that 
was about normal to the flightpath (80° to 90° from the reference axis); and at an 
intermediate angle of about 60°. 


4 





At the outset of the flight program, it was believed that the maneuvers gener- 
ated by the lateral-control input alone would be sufficient to provide most of the 
data for extracting the rotary and sideslip stability derivatives. However, in some 
instances the parameter-estimation calculations resulted in considerable dispersion 
in estimated parameter values. This was apparently largely a result of inadequate 
excitation of the yawing motion. It was concluded that maneuvers generated by a 
rudder input would also be needed — not only for obtaining rudder effectiveness data — 
but also for use with the calculations of the other derivatives. 

A smaller number of maneuvers were performed in which both a rudder and a wheel 
input were employed. In these, the pilot (1) applied a rudder doublet; (2) continued 
with a control-fixed Dutch roll, as in the other tests; and (3) then reestablished 
the initial conditions and applied a wheel doublet. Data were recorded through a 
second Dutch roll. A typical set of time-histories from one of these maneuvers is 
illustrated in figure 6. 


STABILITY-DERIVATIVE PREDICTIONS 


Predicted lateral-directional stability and control characteristics, which 
furnish part of the information used to develop a simulation model for the research 
aircraft, were derived from wind-tunnel tests and theory and from semiempirical rules, 
such as those discussed in references 4 and 5. A large-scale model having an aug- 
mented jet flap and other design features similar to those of the airplane was tested 
in the Ames 40- by 80-Foot Wind Tunnel (refs, 6 and 7). Static stability derivatives 
and lateral-control effectiveness predictions were based on these tests, with appro- 
priate modif ications to account for differences in the areas of the various (control 
and flap) surfaces. 

A method for calculating lateral-directional rotary derivatives from theory for 
a wing with a jet flap is developed in reference 8. An extension of this work by 
J. Farbridge of de Havilland Aircraft of Canada, applicable to the Augmented Jet Flap 
Research Aircraft configuration, is outlined in appendix B. This provides a calcu- 
lated contribution of the wing that, when combined with the estimated contributions 
of the other aerodynamic components, furnished the expressions that were used to 
obtain the rotary derivatives for the simulation model of the complete configuration. 
These expressions and those for the static stability derivatives are given in appen- 
dix C; they are reproduced from reference 9, which describes the simulation model in 
detail, and includes predicted lateral and directional control characteristics. The 
control characteristics are in tabular form and are generally nonlinear. However, 
within the range of rudder deflections used in the tests reported herein, the pre- 
dicted rudder characteristics are very nearly linear with deflection angle. 

Possible effects of the direct hot thrust and the thrust-vector angle on the 
lateral-directional stability derivatives were neglected in these predictions. 


EXTRACTION OF THE STABILITY DERIVATIVES 


The lateral-directional stability and control derivatives were calculated from 
time-history records of the flight maneuvers, using the parameter-estimation procedure 
described later in this section. The application of parameter-estimation methods to 
extract stability derivatives has been discussed and demonstrated in many reports 
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(refs. 10-15). In these procedures, a mathematical model, assumed to represent the 
dynamic behavior of the flight vehicle, generates a response to an input that is a 
measurement of the pilot’s control input. The parameter-estimation calculations 
attempt to determine the set of parameter values for the model that for this input 
causes the model response to best match the measured response. 

The following equations are assumed to represent the lateral-directional dynamics 
of the aircraft in the body-axis system; I x and I z are the aircraft moments of 
inertia about the x and z axes, I X2 is the product of inertia, and a dot over a 
variable denotes a time-derivative. The equations are for a rigid airplane, assuming 


no excitation of the longitudinal dynamics (ref. 13): 

v = wp - ur + g cos 0 sin $ + Ay (1) 

i - *■* ( 2 ) 

«) 

4> = p + r cos <$> tan 0 (4) 


In these expressions, u, v, and w are the components of the true airspeed in the x-, 
y-, and z-direction; p and r are the roll and yaw rates; 0 is the pitch angle; and 
$ is the roll angle. Lateral acceleration Ay and the yawing- and rolling-moment 
terms (Jr, SB') can be represented by the following relations: 

h - z v± + ®° <5> 

< 6 > 

*'IX 5 i + i? o < 7 > 

1 


The terms , and are the dimensional side-force, yawing moment, and 

rolling-moment derivatives, and the terms denote the variables v, r, p, v and 

the control variables defined below. The terms ^ Q , Jf Q9 and are the measurement 

bias terms. With the approximations and substitutions given in appendix D, the state 
equations can be written as follows. 


V 


’ Y v 

Y r - U 

Yp + w 

G 

V 


Y wa 

Y wb 

Y 6 

y r 


^wa 


Y o 

r 


N v 

N r 

N P 

ND 

r 


N wa 

N wb 
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n r 


6 wb 


N o 

P 


L v 

L r 

L P 

LD 

P 

+ 

^wa 

L wb 

L 6 

l r 


6 

+ 

L 

0 



0 

tan 0 

1 

0 

A 


0 

0 

0 

0 


. 6 r. 


0 


( 8 ) 


in which the parameters represent the dimensional derivatives after applying the 
transformations in appendix D, equations (D7) and (DIO) through (D14). The angle 0 
is the constant time-averaged value of the measured pitch attitude during the 
maneuver. Also, u and w, the axial and normal components of velocity, are assumed to 
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to be constant and are given by the following relations in which V av and a av are the 
average speed and angle of attack: 

u = V cos a w = V sin a 

av av av av 

Also in equation (8), G is the product g cos 0, and the parameters ND and LD , 
defined by equations (D12) and (D14) , result from inclusion of v as a variable to 
account for aerodynamic forces at the tail that are produced by sideslip-induced 
sidewash. They were estimated separately (by methods described in ref. 5) and 
entered as fixed values in the equations. 

The displacements of the lateral control surfaces are included as two inputs, 
both defined as effective control-wheel angles. One, 6 wa , is calculated from the 
aileron deflection, and the other, <5^, is the average of wheel angles computed from 
the spoiler and the choke deflections. These latter two sets of surfaces moved 
essentially together but lagged behind the aileron motion. When two inputs were used 
to model the lateral-control displacement instead of a single input representing the 
control-wheel displacement, the computed responses consistently showed better agree- 
ment with measured time-histories. However, the lag of the second input was too 
small to permit the estimation of its separate contribution, and so fixed values were 
entered for one of each of the two lateral input parameters, usually for Y w ^ , N wa , 
and l^b . The values of the lateral-control derivatives shown in the figures in this 
report were obtained by summing the fixed and the extracted values. The third input 
variable, 6, is the time-derivative of S w b , the spoiler-choke displacement, and is 
included to take into account an effect of induced flow at the tail caused by changes 
in circulation resulting from deflection of the lateral control surfaces. The fourth 
input, denotes the rudder displacement. 


Parameter Estimation 

A parameter-estimation program developed by R. Bach (at Northeastern University 
and Ames Research Center) was used to obtain the results presented here. The program 
is described in reference 16, part of which is reproduced in appendix D. It includes 
a provision for selecting a linear-regression or a maximum- likelihood computational 
procedure. The regression solution fits linear combinations of measured states to 
state derivatives and assumes that the state measurements are noise-free. This solu- 
tion required roll and yaw angular acceleration records, p and r, in addition to the 
lateral acceleration. Ay, that was measured. The angular accelerations were calcu- 
lated as the slopes of the curves representing least-squares fits to the roll and yaw 
rate time-histories . 

The maximum-likelihood solution, with state estimation, utilizes a parameter- 
initialization scheme to obtain a first approximation for the response of the airplane 
to the control inputs. Then, in successive steps, the solution iteratively varies 
the parameters to produce a progressively better match of the computed and measured 
responses. Improvement in the matching of the responses in each step is expressed as 
a decrease in an error index, or "cost function," J, defined by equation (D24) in 
appendix D. 

It can be shown that linear-regression methods will cause a bias in parameter 
estimates as a result of measurement noise (see ref. 13). A comparison of results 
obtained with the two procedures, regression and maximum likelihood, offers an indi- 
cation whether such a bias is significant in the present experiments. Evidence 
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presented here indicates that this bias is small. In addition, the linear-regression 
solution yielded parameter estimates that consistently displayed less dispersion than 
those from the maximum- likelihood solution that was employed. For this reason, the 
linear-regression solution was used to obtain most of the results presented in this 
report. Some comments on the nature of the deviations in the estimates produced by 
the maximum-likelihood procedure are included in the discussion of results. 

To process the data, the time-history records of the variables were sampled at a 
selected rate, usually at 50 frames per second (i.e., half the rate at which they were 
recorded) and then, with a filtering routine, the number of data points per variable 
was further reduced to conform to a program array size limitation of 300. A constant 
bias was applied to each measurement to set all of the initial values to zero. The 
parameter-estimation solutions yield estimates of the matrix elements in equation (8) . 
Equations (D15) and (D16) (appendix D) were then used to convert these estimates to 
the nondimensional stability derivatives shown in the figures that present the 
results . 


Sidewash Due to Sideslip 

As mentioned above, the mathematical model for the aircraft motion, equations (1) 
through (7), includes the variable v to account for an effect of induced sidewash 
in the region of the tail. Stability derivatives with respect to this variable were 
not estimated, but as the derivation in appendix D shows, its presence can influence 
the estimated values of the other derivatives. All of the derivative values from 
the tests (but not the predicted values) include incremental contributions due to 
this induced sidewash effect. The magnitude of these increments has been calculated, 
based on predicted values of the £ derivatives given in appendix C (from ref. 9). 
These are listed in appendix E for three sets of flight conditions corresponding to 
angles of attack of 2° (90 knots), 7° (75 knots), and 12° (65 knots). The incremental 
changes are all quite small, except for yaw damping, C n ^. If C nr is corrected by 

the amount shown, the yaw damping is reduced by between 6 and 9%, depending on the 
test conditions and on the actual value of C n ^ within the range of deviation of the 
experimental results. 

RESULTS 


Stability Derivatives 

Figures 7 through 9 present the lateral-directional stability derivatives as 
functions of angle of attack. As indicated in the figure legends, the various symbols 
correspond to different ranges of the jet-momentum coefficient, Cj. The solid symbols 
are results from maneuvers in which the pilot applied a rudder-doublet input followed 
by a control-wheel doublet. Flagged symbols identify results of calculations in 
which it was necessary that two of the rotary derivatives be preset at fixed values. 
The fixed values have been plotted with a smaller symbol size than the computed 
values. Most of these results are from tests produced by rudder doublets, for which 
Lp and N p were fixed. The remaining sets of flagged symbols represent a small frac- 
tion of the tests produced by control-wheel doublets, for which the calculations were 
done with L r and N r set at fixed values. 

Figure 7 shows results for the configuration in which the engine propulsive 
thrust nozzles are at an angle of 8° from the fuselage longitudinal axis; figure 8 is 
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for a 60° angle, and figure 9 is for angles in the range from 80° to 90°. The side- 
force derivatives, Cy^ , Cy^ , and Cy^ (figs. 7(c), 8(c), and 9(c)), are all from 

control-wheel maneuvers. These derivatives from the rudder-excited maneuvers dis- 
played excessive dispersion and are omitted. 

Figures 7 through 9 show, in addition to the experimental results, solid-line 
curves representing (at two values of C j , 0.1 and 0.4) the values for the derivatives 
that were predicted and used in the simulation model, as discussed in the section on 
stability-derivative predictions. The predicted values for and Cy^ are inde- 

pendent of Cj. 


Lateral and Directional Controls 

Rolling-moment, yawing-moment, and side-force derivatives with respect to the 
control-wheel angle are shown in figure 10. As described in the section on reduction 
of data, the control-wheel angle used to define these derivatives is calculated from 
the measured displacements of the lateral-control surfaces, with gearing ratios from 
ground calibrations. Also shown in the same figures are the derivatives with respect 
to the rate of displacement of the lateral control. As indicated in the discussion 
of the mathematical model for the aircraft dynamics, these derivatives are included 
to account for forces and moments that result from downwash and sidewash (at the tail) 
induced by the lateral-control deflection. 

In figure 11, rolling-moment coefficients are shown as a function of control- 
wheel angle. The two solid-line curves are predictions used in the simulation model, 
based on large-scale wind-tunnel tests, for two values of C j , 0.1 and 0.35. The 
cross-hatched region corresponds to the results from the parameter-estimation program, 
which yields the slopes of linear curves through the origin. The region shown was 
obtained by averaging the results of all tests for which the average angle of attack 
was in the range from 2° to 10°, its width being computed as twice the standard 
deviation. 

Figure 12 shows the directional control characteristics. These are the deriva- 
tives with respect to the rudder-deflection angle measured (in degrees) on the forward 
panel of the two-panel rudder. The rudder characteristics used in the simulation 
model are shown as the solid- line curves and apply to a range of rudder-deflection 
angles within which the predicted variation is linear, that is, up to about 15°. 

For the flight experiments, the deflections were within this range. 


DISCUSSION 


The test results (figs. 7-9) indicate that the magnitude of the roll-damping 
derivative (C^ ) is significantly greater than that used in the simulation model and 

that at angles of attack near zero, it varies more rapidly with angle of attack than 
predicted. The large estimated roll damping indicates that the section lift-curve 
slopes for the outer portions of the wing are relatively high. The usual deleterious 
effects of large aileron and spoiler deflection on lift-curve slopes appear to be 
alleviated by boundary- layer blowing at the aileron and by the leading-edge slat 
configuration that is used. 
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The graphs of show a large dihedral effect at small and negative angles of 

attack and a gradual decrease in this effect as the angle of attack becomes large. 

The curves representing the predicted are based on wind-tunnel measurements and 

indicate values near zero for all conditions. The wind-tunnel model had no geometric 
dihedral, whereas the flight vehicle was constructed with a 5° dihedral in the wing 
panels outboard of the engine nacelles. The graphs showing C n g indicate a decrease 

in this derivative with angle of attack, compared with a predicted constant value, 
but the levels at moderate angles of attack are in agreement. The levels of yaw 
damping and the derivative generally agree with predictions; C n ^ is in agree- 

ment at small angles of attack but shows a somewhat greater rate of variation with 
angle of attack. 

Side-force data plotted in figures 7(c), 8(c), and 9(c) indicate that the values 
of derivatives Cy^ and Cyg are near the predicted levels; there is little variation 

of Cyg with angle of attack, although a negative variation was predicted. The 

average of the estimates for Cy r is larger than predicted by 30% to 40% at moderate 

angles of attack and by over 80% at angles above 8°. 

Among all of the derivatives shown in figures 7-9, little or no consistent effect 
of variation of jet-momentum coefficient Cj is evident. Also for all the deriva- 
tives, there is significantly more dispersion in the estimates from tests in which 
the engine thrust nozzles were at angles in the range from 80° to 90° than for the 
smaller angles (8° and 60°). At angles of attack where comparisons can be made, the 
effects on the derivatives of changing the angles at which the nozzles were set 
appear to be small. 

Figure 10(a), which shows the lateral-control effectiveness Cj^ for the con- 
figuration with the engine thrust nozzles at 8°, indicates an increase in effective- 
ness with angle of attack. With the nozzles at 60° and 80° to 90°, a consistent 
effect of angle of attack is not apparent. For all of the nozzle angles, the roll- 
effectiveness data show little consistent effect of variation of Cj. 

From the graphs showing the lateral-control rate derivatives (fig. 10) , it is 
seen that the estimated values of the yawing-moment derivative C n ^ are predominently 

negative in sign. This can be interpreted as indicating that a positive lateral- 
control deflection induces a sidewash at the vertical tail, causing an added, but 
delayed, positive yawing moment. Estimated values for the rolling-moment derivative 
C££ are also mostly negative in sign. This means that part of the rolling-moment 

response to the lateral control is delayed, an effect that can be attributed to an 
induced downwash that varies along the span of the horizontal tail, with its strength 
increasing toward the right for a positive control deflection. 


Rolling Moment versus Wheel Angle 

In figure 11 (rolling-moment coefficient plotted as a function of control-wheel 
angle) , the data from parameter-estimation calculations (the cross-hatched region) 
represent results from 38 maneuvers covering a range of Cj from 0.1 to 0.38. As 
mentioned earlier, there was no consistent correlation of (the slopes of the 
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curves represented by this region) with this Cj variation. This differs from the 
predicted effect of C j , illustrated by the solid line curves in figure 11, where 
there is a large increase in effectiveness with an increase in Cj from 0.1 to 0.35. 
The results from the flight experiments indicate a considerably greater rolling 
moment due to control deflection than that predicted for the smaller Cj. The flight 
results also show larger rolling effectiveness than predicted for Cj of 0.35 at the 
larger control-wheel angles. The broken-line curve shown in figure 11 represents 
rolling-moment coefficients calculated from data reported in reference 3; the data 
were obtained from roll-reversal maneuvers, in which the roll acceleration was deter- 
mined at times when the roll rates were near zero. These results indicate a somewhat 
larger roll effectiveness than that from the parameter-estimation calculations for 
control-wheel angles in the range from 20° to 40°; at smaller and larger wheel angles, 
however, the results from the two techniques are generally in agreement. The 
parameter-estimation results are from maneuvers that were generated by applying a wide 
range of control-wheel time-history functions. The maximum control-wheel angles used 
for these input functions ranged from 13° to 62°. For approximately half of the 
maneuvers, the maximum wheel angles were greater than 40°, and the average of the 
maximum wheel angles for all of the tests was 37°. 


Rudder Effectiveness 

As determined from the flight experiments, the yawing moment produced by the 
rudder (fig. 12) agrees well with that used in the simulation model. The average 
values for the rolling-moment and side-force derivatives from these tests are about 
30% smaller than those in the simulation model. In reference 9, a relation between 
C n ^ R and Cy^ R is obtained by assuming that the rudder yawing moment is equal to the 

product of the side force and the distance from the aircraft center of gravity to the 
aerodynamic center of the vertical tail. The data from the flight experiments indi- 
cate that the rudder contributes a significant additional yawing moment attributable 
to an aerodynamic couple on the vertical tail. 


Maximum-Likelihood Results 

In figure 13, stability-derivative estimates obtained using the maximum- 
likelihood solution (appendix D) are compared with estimates from the linear- 
regression solution. These results are all from flight maneuvers in which the pilot 
first applied a rudder doublet and then a control— wheel doublet. Data represented by 
the plain symbols are those shown earlier (with solid symbols in figs. 7-9) from 
linear regression, and the flagged symbols are from the maximum- likelihood computa- 
tion. As before, the various symbols correspond to different Cj ranges. A notable 
difference in the results obtained with the two computational procedures is the 
increased dispersion in the maximum— likelihood data. The increase is moderate in the 
results for the aircraft configuration in which the vectored thrust nozzles were at 
8° and 60° (fig. 13(a)) and is quite large for the 85° nozzle angle (fig. 13(b)). 

It is generally recognized (and shown, for example in ref. 17) that there will 
be biased estimates resulting from linear-regression solutions, if measurement noise 
has an important effect. Some indication of such biases can be discerned from 
figure 13, although scatter in the data prevents an accurate evaluation of their 
magnitudes. The regression solutions yield slightly larger values for C n ^ and 

indicate slightly less roll damping than the maximum- likelihood solutions, but these 
differences are quite small. Because of this small effect and because using linear 
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regression not only resulted in significantly less scatter but also entailed consid- 
erably less computer execution time, it is clearly the preferred method for process- 
ing the data from the tests reported here. 1 Situations exist, however, in which the 
maximum- likelihood method with state estimation might be chosen instead; for example, 
when measurements of all the states are not available, or when measurement noise is 
so large that it causes significantly biased estimates. For the maximum- likelihood 
procedure to offer a suitable alternative, it is important that the large dispersion 
be reduced. A reduction in this scatter was achieved by some minor changes in the 
way the computations were handled, as discussed below. 

The comments that follow are based on the results of applying the maximum- 
likelihood procedure to all of the relatively large number of time-history records 
that are the subject of this report. The largest deviations in the parameter esti- 
mates were not random in nature. It was evident that when these deviations occurred, 
an apparent error in one parameter was related to errors in one or more other 
parameters. In some cases, such correlated errors can be traced to the fact that 
during the latter part of the maneuvers, when the motion of the airplane is a control- 
fixed Dutch roll, the side velocity, roll rate, and yaw rate are linearly related. 

For such motion, errors in one of the three derivatives C£ , r , and Cj^ (or C n ^ , 

C nr , and C n g) are compensated for by errors in the other two, with the result that 

the solution will not necessarily converge to the correct parameter values. Presum- 
ably, because of the presence of random perturbations of the aircraft's motion caused 
by turbulence, the earlier portion of the maneuver does not contribute strongly enough 
to assure a unique solution. By fixing the value of one of the three derivatives at 
its correct value, these errors in the other derivatives are minimized. When errors 
of this type were evident in the results from maneuvers generated by a rudder doublet 
followed by a control-wheel doublet, better results were obtained by dividing the 
data records into two time intervals and treating them as separate maneuvers. A 
standard procedure was adopted of inserting fixed values for Lp and Np in the solu- 
tions for the rudder generated maneuvers and fixed values for L r and N r for lateral- 
control maneuvers. By alternating the execution of the two solutions and successively 
improving the selection of the fixed values, a consistent set of derivatives could 
usually be determined in a few iterations. 

Maneuvers in which only the rudder or only the control wheel was used to excite 
the motion also were routinely treated in this way when the maximum-likelihood proce- 
dure was used. Sometimes the choice of the values to be fixed could only be made by 
an interpolation or an extrapolation from results at other test conditions. When the 
linear-regression solution was used, and if the maneuver was produced solely by a 
rudder input, the same process of fixing two derivatives (Lp and Np) was necessary; 
but in most instances when maneuvers were generated by a control-wheel input, the 
calculations yielded satisfactory results without fixing any of these four derivative 
values. These data are indicated in figures 7-9 by open symbols without flags; the 
data for solutions requiring fixed values are identified by the flagged symbols. From 


Reference 14 presents the results of applying a linear-regression computation 
to obtain longitudinal stability derivatives for the airplane that is the subject of 
this report. These results are compared with results from a quasi-linearization 
(modified Newton-Raphson) method, which, like maximum likelihood, is an iterative, 
"output-error" procedure, with state estimation. This reference notes that the 
regression method provided smaller standard deviations in the estimated parameter 
values and better agreement with wind-tunnel measurements than did the quasi- 
linearization method. 
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the preceding comments it is seen that when the maximum- likelihood solution was used, 
for some test conditions it was necessary to have available flight records from a 
greater number of maneuvers (i.e., additional rudder-excited maneuvers) in order to 
obtain a satisfactory set of rotary and sideslip derivatives than would be required 
by the linear-regression solution. 


A significant fraction of the maneuvers for which maximum- likelihood solutions 
were obtained produced results which showed evidence of "nonrandom" deviations of 
another type. 2 These occur when an error in one of the roll-rate derivatives 
(C^, C n ^, or Cy^) is related to an error in the corresponding lateral-control deriv- 
ative (Cj^, c n w » or Cy w ) . The size of this type of error was observed to be adversely 

influenced by the following factors associated with the test conditions: (1) a high 

level of engine thrust at low airspeed, so that the thrust coefficient and Cj were 
large; (2) a setting of the engine thrust nozzles at angles in the range of 80° to 
90°; and (3) a diminished angular range through which the lateral-control surfaces 
were moved. With regard to the last of these, there were significant deviations in 
the estimates from a majority of the solutions for flight data in which the maximum 
displacement of the control wheel was less than 25°. 


An example of such errors is observed when, in the same solution, both the cal- 
culated value for the lateral-control effectiveness Cj^ and the calculated value for 

roll damping are too large. When solutions showed evidence of such errors, the calcu- 
lations were repeated with fixed values for the lateral-control derivatives. Because 
the flight experiments covered an extensive range of test conditions and because most 
of the solutions were not subject to these errors, there was generally sufficient 
information to permit the selection of suitable values for the lateral-control 
derivatives . 


Results of applying the maximum-likelihood solution to the data from the maneu- 
vers produced wholly or in part by a lateral-control input when the thrust nozzles 
were set (nominally) at 85°, are plotted in figures 14(a) and 14(b). On the left in 
each figure are results from solutions in which the three lateral-control derivatives, 
C^, C n ^, and Cy^, were among the set of estimated parameters. On the right, results 

from the same test data are plotted but with the values for one or both of the two 
control derivatives, and C n ^ 9 fixed, if in the earlier solution they differed 


substantially from the mean of the results from the other tests. The flagged symbols 
in these figures identify the lateral-control derivatives that were set at fixed 
values. 3 The amount of scatter in the graphs of and C n was significantly 

reduced by inserting the fixed quantities. P ^ 


Reference 18 presents and discusses other examples illustrating correlation of 
input and response parameters. In those examples, the correlation is a result of the 
relationship of elevator displacement and the pitch rate. As in the case of the 
present investigation, results obtained using a maximum- likelihood algorithm in the 


2 Use of the linear-regression solution did not eliminate deviations of this type, 
but it did reduce the number of instances in which they appeared and decreased the 
magnitude of the deviations. 

3 In figure 14, solid symbols represent maneuvers that included rudder inputs; 
open symbols indicate maneuvers made with control-wheel inputs. The latter solutions 
were obtained with L r and N r set at fixed values, selected using the derivative 
values shown in figures 9(a) and 9(b) as a guide. 
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examples of reference 18 exhibit larger deviations attributable to this type of corre- 
lation than those obtained with linear regression. 

CONCLUDING REMARKS 


A linear-regression parameter-estimation procedure was used with data from flight 
experiments to obtain lateral-directional stability and control derivatives for the 
Augmented Jet Flap Research Aircraft. The tests, which were conducted with the jet 
flaps deflected to a position typical of a landing approach (65°), covered a large 
range of angles of attack below the stall, as well as a range of engine power levels 
up to the rated maximum for continuous operation. The aircraft design enables the 
pilot to control the direction of the engine propulsive thrust by moving the jet- 
exhaust nozzles through an angular range of over 90°. The tests included a study of 
the effects on the derivatives of setting the nozzles — in addition to their positions 
for forward thrust — at a downward angle that was approximately normal to the aircraft 
longitudinal axis, and at an intermediate angle that was 60° from the axis. In 
general, variations in the angle at which the nozzles were set and in the flap jet- 
momentum coefficient had a rather small effect on the levels of the stability- 
derivative estimates. 

The yaw damping (C n ^) and the derivative C^ from the flight experiments are 

generally in agreement with predictions that had been made before the flight program. 
The flight results also indicate that there was significantly more roll damping at 
all angles of attack than had been predicted, and that at small angles of attack the 
effect on roll damping of varying the angle of attack was greater than predicted. 

The derivative C n ^ also shows a variation with angle of attack that is somewhat 

greater than predicted. Values for the directional derivative C n ^ are in agreement 

with the predicted levels at moderate angles of attack, but they display a sizable 
decrease with angle of attack that was not predicted. 

The results obtained using a maximum-likelihood parameter-estimation procedure 
are compared with results from the linear-regression calculations. For test condi- 
tions corresponding to large values of the flap jet-momentum coefficient (Cj) , 
reduced airspeeds, and to a setting of the engine vectored-thrust nozzles at large 
angular positions (between 80° and 90°) , the maximum-likelihood solutions produced 
significantly greater dispersion in the estimated-parameter values than did the linear 
regression. The comparison also indicated that if the linear-regression estimates 
are biased, as a result of the effects of measurement noise, for these tests the 
amount of such bias is not large. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, August 24, 1982 
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APPENDIX A 




SIDESLIP VANE ANGLE CORRECTIONS 


The sideslip angles indicated by a vane on the nose boom were corrected for, 
sidewash at the vane, using the following relations and assuming that the error 
varies linearly with the sine of the sideslip angle. Defining as the indicated 

sideslip angle transferred to the airplane center of gravity and 3 C as the actual 
sideslip angle, a correction factor Ky is defined as 

K v = sin 3 c /sin g ± (Al) 

If K v is known, a corrected side velocity v c can be obtained with the relation 

v = K V sin £. 
c v 1 

The longitudinal and vertical components of the true airspeed V, in the body-axis 
system (at angle of attack a) are, respectively, 

u - V cos 3 C cos a 


and 


w - V cos B c sin a 

The rate of change of side velocity is given by the equation 

dv /dt = A - ur + wp + g cos 0 sin <j> (A2) 

c y 

The terms Ay, r, p, g, 0, and cf> are defined in the symbol list. Values of K v for 
various flight conditions were determined by performing least-squares fits of the 
rates of change of the side velocities computed from the vane measurements to the 
time-histories of dv c /dt calculated from equation (A2) . 
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APPENDIX B 


I IHIIIIIIIIIIIIIIIIIIIII 


ROTARY-DERIVATIVE PREDICTIONS 


This appendix outlines work by J. Farbridge (of de Havilland of Canada) to modify 
and extend the results of reference 8, which gives expressions for calculating the 
lateral-directional rotary-stability derivatives for wings with jet flaps. The 
derivations in reference 8 include the assumptions that the chord of the flap is 
small (approaching zero) and that the flap deflection angle is small. Farbridge 
incorporates the results of reference 19 to account for the effect on the derivatives 
of a larger flap chord and uses results from reference 20 to estimate the added effect 
of a large flap deflection. In reference 20 the effect of the flap angle on the 
lift coefficient of an airfoil section without flap blowing is calculated from 
theory. By analyzing and correlating experimental lift data on wings with jet flaps, 
Farbridge deduced the following relation for the lift coefficient of an airfoil due 
to jet flap blowing: 

C Lco = 3.5(0,, + 6 f ) - 1.08(0* + R°- 25 6 f )c;- 75 (Bl) 

where R c is the ratio of flap chord to Wing chord and a w is the wing angle of 
attack. The total section lift coefficient Ct is defined as the sum of Ct t and 

Cl > the lift on a wing without blowing. The lift coefficient for a finite span 

°°o 

wing of aspect ratio A is defined by the following relation: 

C L = (ttA + 2C j )C l /B (B2) 

A 


where 


B = ttA + BCl^/Soi + [ C j / ( 6 ^ + a w ) ] BCl^/BCj 


(B2) 


The following expressions illustrate the changes in the equation for roll damping 
that result from Farbridge T s modifications. From reference 8, 


r _ _L / jrA , 4 \ 9Cl °° [ 

c «. p 2 V 4 + 3 C jy 9a [j 


3C t 


ttA + 2 


3a 




(B3) 


In this equation, is replaced by the relation 


C L = C L + C L 

^oo ^coj -q 


in which Cl^j is from equation (Bl), and is from airfoil-section theory and 

reference 20. This, with some rearrangement, becomes the following: 


ttA 

8 


2.42 C 0 / 75 + 

J 


1 8 “ v„ 


1 + (16/3) (Cj/ttA) 


(B4) 


16 


min ii iiMM i ■iiiii 


i ■!!■■■ min 


where 


F = ttA + 2 




(1 - 0.3086 R°* 2S )6 f + 0.690 a + 1.726 


84 + 5.25 


6^ + a + 2. 5 


and 


a = a - 2.5 
w 

When the expression for C n ^ from reference 8 is similarly combined with equa- 
tion (Bl) (and its partial derivatives), the following relation is obtained: 


n 


8 L a ttA 
A 


TW_ ' 1 - 3^2.42 C*-» + (S) Cj j - 

lb 


+ (32/9) (Cj/ttA) 


) 


(B5) 


where Cl is the lift coefficient of the wing of the flight vehicle. The expression 

A 

for from reference 8 with these modifications is 


= C 


ttA 


([1 + (16/3)C J /irA]K 


L a 8 (ttA + 2Cj) 


+ 1 


(B6) 


where 




K = hrA + 2.42 C 


/3C t \ 57.3 C T 1 

• 7 5 + ( Z»l + . ^ 

\ 3a / C 2 ( 6 f + a + 2.5) 


(1 - H) + I 


H = 5.25 C°- 75 


(1 - 0.3086 R ° ‘ 2 5 ) 6 f + 2.416 a + 6.04 


I = 3.624 C°* 75 + 


HC t + 0.405 C° t - 75 [(R°- 25 - 3.241)6 f - 2.237 a- 5.592] 
^oo J C L 


6^ 4- a + 2 . 5 


The contribution of the wing to the damping in yaw Cn^, as indicated in refer- 
ence 8, is computed from wing-section drag and is separated into two parts: one due 

to profile drag C n and one due to induced drag Cn . • As in reference 8, 

r 0 r i 

c n rQ = ~ (1 / 4 )Cd 0 > where Cq o is the profile-drag coefficient. The modified expres- 
sion for the contribution due to induced drag is 


0 ttA(1 + 32 C / 9ttA) __ 

P ^ P 2 J 

° n r . - 8 a ~ N *, F 

i A (ttA + 2C ) 2 


(B7) 


The derivatives given by the above equations are referred to stability axes. 



APPENDIX C 


STABILITY DERIVATIVES USED IN THE SIMULATION MODEL 

The following formulas for lateral-directional stability derivatives are 
from reference 9. They are part of a computer program for the dynamic simulation of 
the Augmented Jet Research Aircraft. In these expressions, the derivatives are 
referenced to stability axes, whereas those in the figures (the "predicted” deriva- 
tives), like the experimental data, are in the body-axis system. As given here, they 
apply only for a 65° deflection of the landing flaps. The angle of attack in the 
formulas is measured from the aircraft fuselage reference and is expressed in degrees. 

C = -1.166 - 0.0344 a 

y e 

C = -0.0157(sin a + 4.037 cos a) (Cl) 

y h 

C = 0. 1828(sin a + 4.037 cos a) - 0.0112 C 2 
y r l A 

C = 0.216 C 0 - 0.15 C - 0.1828(cos a - 4.037 sin a) + 0.186 
y ^ l A 

P P A 

C = 0.000766 a - 0.0331 C T + 0.251 

n e J 

C = 0. 00193(sin a + 4.037 cos a ) 2 (C2) 

n 6 

C = 0.0735 - 0.0119 C 2 (1 - 0.649 C 2/2 )(l + 0.00275 a) 
n L . J 

r A 

- 0.0225(sin a + 4.037 cos a ) 2 - 0.0584 

C = -0.00578 - C T [0.0860 + 0.00053 a - 0.02 C T ] 
n L 3 

P A 

+ 0.0225(sin a + 4.037 cos a) (cos a - 4.037 sin a - 1.021) 

C Q =0 (from tests of wind tunnel model with no dihedral) 

0 

C = -0.00193(sin a + 4.037 cos a) (cos a - 4.037 sin a) (C3) 

C = 0.239 C T (1 - 0.285 C 2/2 )(l + 0.00231 a) 
r L A J 

+ 0.0225(sin a + 4.037 cos a) (cos a - 4.037 sin a) + 0.045 

C £ = 0.00440 a - 0.19 Cj - 0.405 
P 
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APPENDIX D 




MATHEMATICAL MODEL AND PARAMETER-ESTIMATION EQUATIONS 


This appendix presents the relations and assumptions used to obtain equation (8) 
and describes the parameter-estimation procedure used in calculating the stability- 
derivative estimates in this report. 


Derivation of the State Equations 


The summation terms in equations (5)- (7) represent the linearized aerodynamic 
force and moment derivatives, as follows. 


E l?. = ^v + ^v+^r + ^p+^ 6 + ^ u + 

i \r v r ^p^ ^wa wa wb wb o R R 

V' JV E . = jy v + jV-v + jV r -V jV p + jY 6 + JV ,6 , + ^6 + 

/ * T v v r p r wa wa wb wb 6 R R 

£ . = g v + + &r + c?p+g> 6 + i? , 6 , + ^6 + 3? 6 

i v v r p r wa wa wb wb o R R 

With the following definitions 


(Dl) 

(D2) 

(D3) 


A = I /I 
xz z 


B = I /I 

XZ X 


(D4) 


C = 1 - I z /(I I ) 

xz x z y 


eliminating p and r, respectively, from equations (2) and (3) yields the expressions: 


rC = v(^ + A5? v ) + v(^ + A^) + r(j^ + A^ r ) + pOt + A5? p ) + + ASp 


i = i 
4 


(D5) 


pC = v(5? v + B«/T) + vCS?^ + &>) + r(5? r + B^.) + p (SC p + B*/^) + ^ & ±^± + Bt ^. ) 


where 6 j_ denotes the control variables, 6 wa , 6 W ^ , 6, and 6 r. From equation (1), and 
defining D = 1 - 


v = a 


v = | + (®^ - u) r + 4- w)p + g cos 0 sin (j) + 


Z Vi / D 

i=i x J/ 


(D6) 
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Then with the substitutions 


N 5 = ^ + A ^ ? )/c 

and 

= 0 ? + M)/C 

and v from equation (D6) , equation set (D5) can be written as 

r = v (N^ + N^ v /D) + r[lT. + N^(» r - u)/D] + p[N^ + N£(* + w)/D] 

4 

+ 53 6 i( N s + ^ + cos 6 sin <f) ^ D 

i=i i i 

p = v(L; + L^ v /D) + r[I/ + L^(» r - u)/D] + p[I/ + L^(«r + w)/D] 

4 

+ 53 6 i( L 6 + /D) + ( L ^g cos 6 sin <f>)/D 

i=i i i 


(D7) 


(D8) 


(D9) 


Because roll angles during the tests remained relatively small, for this analysis it 
can be assumed that sin (J) = (p and that <j> = p + r tan 0 a , where an average (con- 
stant) value, 0 a , is used in place of the varying pitch angle. The state equations 
can then be written as in equation (8) with the parameters defined by the relations 
below. 


Y 

v 

Y + w 
P 


N 

v 

N 

P 

ND 

L 

v 

L 

P 

LD 


V D ’ 

(» p + w) /D , 


N' + N ! & /D , 

V V V 


N T + N f . + w) /D , 
P v p 

(N^g cos 0)/D 

l; + L! v° , 

+ L -^p + W >/ D * 

(L^g cos 6)/D 



(» r ~ u)/D , 

*s /D 

i 

N; + N^(^ r - u)/D 
N' + N!^. /D 

O . VO. 

1 1 


(DIO) 

(Dll) 

(D12) 


l; + U(* r - u)/D , 

(D13) 

L 6 . + L ^6 /D 

x i 

(D14) 


A value for D based on Cy£ given by equation (Cl) is estimated in appendix E. 

For a typical set of test conditions, D is calculated to be 1.002, and it varies only 
slightly from this at other conditions of the tests. 
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The following relations were used to compute the stability and control deriva- 
tives that are presented in this report. The true velocity V and the dynamic pres- 
sure q [=(1/2) pV 2 ] are the average values for the time of the maneuver. 

C = Y v VW/qSg) 

p 

C = (N - AL )VI / (qSb) 
n^ v v y z 7 

= (L v - BN v )VI x /(qSb) 

B 

C = 2Y r VW/(qSbg) 
y r 

C n = 2(N r " AL r )VI z /(qSb 2 ) (D15) 

r 

= 2(L r - BN r )VI x /(qSb 2 ) 
r 

C = 2Y VW / (qSbg) 

y p p 

C = 2(N - AL )VI / (qSb 2 ) 

n p p z' 

P 

C £ = 2 (L - BN )VI x /(qSb 2 ) 

Ti ir tr 


The cont rol-def lection derivatives are computed from the following equations, where 
6 represents the variabl es 6 wa , <5 w b > <5 » and Sr. 


c = Y 6 W/(qSg) 

* 6 

C n = (N 6 - AL 6 )I z /(qSb) (D16) 

6 

C . - ' BN s )I x /(5sb) 

0 


The control-wheel derivatives are calculated from the relation Cj\^ = 
which A corresponds to the subscripts y, n, and 


Parameter Estimation 

It is assumed that the system is described by the vector differential equation 

x(t) = Fx(t ) + Gu (t ) + w ( t ) ; x(t Q ) = x Q (D17) 


where x(t) is a state vector, u(t) is an input vector, w(t) is the state (or process) 
noise vector, x Q is the initial condition vector, and F and G are, respectively, 
the stability and control matrices. The system output can be expressed as 


y T (t) 




(D18) 
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where x m (t) and a m (t) are state and state-derivative vector elements; the complete 
state-derivative vector can be expressed as 

a(t) = F x(t) + Gu(t) + w(t) (D19) 

a 

where F a is a subset of F. The parameter-estimation procedure outlined here 
(excerpted from ref. 16) is formulated to allow the selection of a maximum- likelihood 
or a linear-regression solution to estimate the unknown parameters in the stability 
and control matrices. The data system yields a time-history record on an interval 
(t Q »tf) with samples spaced T seconds apart comprising control inputs u(t) and the 
record 


z(t ± ) = y(t i ) + v(t ± ) , i = 0, 1 , . . . , N - 1 

where and v(tt) are measurement and measurement-noise vectors, respectively. 

(Other assumptions: that the measurement noise and the process noise w(t) are white- 

Gaussian with means, respectively, of yb and wb . ) 

In the maximum-likelihood parameter-estimation procedure, a model estimate z(t) 
is matched to a data record z(t) by adjusting model parameters with an iterative 
strategy to minimize an error criterion. If the flight maneuver is performed under 
conditions that are relatively free of atmospheric turbulence, the random component 
of w(t) in equation (D17) can be ignored and a state estimate is obtained from the 
model 

x ( t ) = Fx(t ) + Gu (t ) + w b ; x(t Q ) = x Q (D20) 

where Wb can account for bias error in the measurement of u(t) . The model output 
is given by 


y T( t } = [x*,£*(t)] 

where the state-derivative estimates would be taken as a subset of 

a(t) = F x(t) + Gu(t) + w, (D21) 

a d 


Then an estimate of the data record is formed as 

z(t i ) = y(t i ) + v b ; i = 0, 1 , . . . , N - 1 


(D22) 


Here Vb represents a constant bias in the measurement of y(t). 
The parameter vector can be written as 


T 

P 


T T 
P FG’ P w, 




(D23) 


For the maximum-likelihood formulation the criterion for matching the model output to 
the data record is the "cost function": 
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J = (1/2) 


(D24) 


N-i 

E 

i=o 


E [z(t ± ) - S(t.)]V 1 [z(t.) - z(t i )] + (N/2)log|] 


where the minimizing value of B is given by 


N-i 


B = (1/N) E [zCV = 2(t.)][z(t.) - z(t.)]' 


(D25) 


1=0 


This criterion is minimized with a modified Newton-Raphson (or quasi-linearization) 
algorithm (ref. 11). A parameter change is computed as 6p = -M 1 (3J/3p) T , where M 
is given by 


N-i 


M = £ S i (t.)B 1 S(t i ) 


(D26) 


1=0 


and the gradient vector 3J/3p is 


N-i 


(9J/9p) = -E tz(t.) - z(t i )]V 1 S(t i ) 
i = o 

In the above equations, S(t) is the sensitivity matrix 

S(t) = 3z (t ) / 3p 

For a linear-regression fit with no process noise, 

z(t) = a m (t) 

and the kth column of S(t) can be obtained as a subset of the vector 
[3a(t)/3p k ] = (3F a /Bp k )x(t) + (3G/3p k )u(t) + (3w b /3p k ) 


(D27 ) 


if, as is the case in the present application, all the states are available among the 
measurements. For a maximum-likelihood solution, the evaluation of S(t) requires the 
solution of as many as n e differential equations during each iteration, where n e 
is the product of the number of states and the number of unknown parameters. Thus, 
in the process-noise-free case, from equation (D22) , 

S (t ) = [3y(t)/3p] + (3v b /3p) (D28) 

which implies that both 3x m (t)/9p and 3a m (t)/3p must be computed. The kth column 
of the state-sensitivity matrix 3x(t)/3p is obtained from equation (D20) as the 
solution of the vector differential equation 

[9x(t)/9p k ] = F[9x(t)/9p k ] + (9F/9p k )x(t) + (9G/9p k )u(t) + 9w b /9p k (D29) 


with initial conditions 
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0 


3x(0)/9p k 


, p, not in x 
r k i 

e , p, in x 

r lr 


Here e is a vector with zero elements, except for a one in the row corresponding to 
an initial-condition parameter p^. 

To increase the convergence range of the quasi-linearization algorithm, a pro- 
cedure known as a "Denery start" (ref. 17) is included in the computational program. 
In this procedure, equation (D29) is modified by solving the state sensitivity equa- 
tions, using state measurements in place of x(t). 
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APPENDIX E 


COMPUTED EFFECTS OF B DERIVATIVES ON PARAMETER ESTIMATES 


In this appendix the effect on the stability derivatives of sidewash due to 
sideslip is calculated. The extracted parameter values include contributions from 
this induced sidewash. These contributions are shown in equations (DIO) through (D14) 
as the terms in which the subscript v (or the symbol D or both) appears. In the 
expressions below, the magnitude of their contribution to the stability and control 
derivative values is computed, based on the predicted values for Cyg, C n g, and 

that were used in the simulation model and are presented in appendix C. To apply a 
correction to the extracted derivative values for this effect, the following incre- 
ments would be subtracted from values shown in this report. 

Side force: 


AC -EC C 


y X y g 


(El) 


where E = pSbg/(4W) and X denotes the variables (3, r, p, 6 W , and 6 r. 
Yawing moment: 


Rolling moment: 


AC = EC C 

n B 


y e n B 


= -[cos a - EC JC 

n r V y r/ n B 

= (sin a + EC )C 

n \ yin* 

P \ p/ 3 


AC 


AC = EC C 


n . 


y 6 n § 


AC £ 

3 

= EC C 
y 6 ‘6 

ac sl 

= - (cos a 

r 

\ 

AC * 

= (sin a • 

p 

V 

AC x 

= EC C. 


+ EC 




The side force increments can all be expressed as a fraction of the side-force 
derivative. 



With Cy^ 


given by equation (Cl), this relation indicates a change of less than 0.4% 


of the derivative values under any conditions of the tests. Also using equation (Cl), 
the value for D in equations (D-6) through (D14) can be computed: 


D = 1.0 - 

v 

= 1.0 - EC (E3) 

y 6 

= 1.002 


for a typical set of test conditions and varies only slightly for other test 
conditions . 

Increments to the yawing- and rolling-moment derivatives from the above formulas 
are listed below for three angles of attack, 2°, 7°, and 12° (corresponding to, 
respectively, assumed equivalent airspeeds of 90, 75, and 65 knots), an aircraft gross 
weight of 178,000 N (40,000 lb) and an altitude of 925 m (3,000 ft). 


a 

2° 

7° 

12° 

AC 

n e 

-0.0014 

-0.0014 

-0.0014 

AC * 

p 

.0003 

.0002 

.0000 

AC 

n 

r 

-.0308 

-.0299 

-.0283 

AC * 

r 

.0066 

.0038 

.0009 

AC 

n 

P 

.0007 

.0033 

.0061 

AC £ 

P 

-.0001 

-.0004 

-.0002 

AC 

n 

w 

.0008 

.0005 

.0004 

AC £ 

w 

.0002 

.0001 

.0000 

AC 

\ 

.0000 

.0000 

.0000 


.0000 

.0000 

.0000 


Except for yaw damping 

5 ^ n r * all of these values represent negligibly small 

corrections to the extracted derivative values. If C n is corrected by the amount 
shown above, the yaw damping would be reduced by between 6% and 9%. 
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1 


TABLE 1.- AUGMENTED JET FLAP RESEARCH AIRCRAFT: WEIGHTS AND AREAS 


Weights, N (lb) 

Maximum gross 213,000 (48,000) 

Operational, empty 145,000 (32,600) 


Areas, m 2 (ft 2 ) 
Wing 

Horizontal tail 
Vertical tail 


80.36 (865) 
21.65 (233) 
14.12 (152) 


Other geometric data 

Wing dihedral, starting 
at station 5.36 m from 
plane of symmetry, deg 
Ailerons 

Span, m (ft) 

Chord aft of hinge 
line, m (ft) 

Spoilers 

Span, m (ft) 

Chord, m (ft) 

Position of hinge line 
% wing chord, 
average 
Flaps 

Span (each side) , 
m (ft) 

Chord aft of hinge 
line, m (ft) 


3.51(11.50) 

.61 (2.01) 

3.44(11.30) 
.36 (1.18) 

62.4 

7.01(23.0) 
.98 (3.2) 


Moments of inertia, kg-m 2 (slug-ft 2 ) 
(Gross weight 178,000 N (40,000 lb)) 
With original elevator control 
(spring tab): 

I x = 357,000 (263,300) 

I y = 278,000 (205,000) 

I z = 587,000 (432,900) 

With modified elevator control 
(powered elevator) 

I x = 361,000 (266,300) 

I y = 316,000 (233,100) 

I z = 620,000 (457,300) 

Product of inertia 

I xz = 48,130 (35,500) 
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[ 4 — 9.3 m 
CONICAL' (30.5 ft) 

NOZZLES 


X ROLLS-RO> 
SPEY MK 8( 


Figure 2.- Three-view drawing of 








L.E. SLAT 




Figure 3.- Details of the engine installation, lateral control surfaces, 

and augmented jet-flap components . 
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FILLED SYMBOL - BOTH RUDDER AND LATERAL 
CONTROL INPUTS. 

FLAGGED SYMBOL - TWO PARAMETERS FIXED. 




(a) C £ , C £ and vs a. 

p r e 

Figure 7.- Lateral-directional stability derivatives: nozzles at 8°. 
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CONTROL INPUTS. 

FLAGGED SYMBOL - TWO PARAMETERS FIXED. 





(a) C , C , and C vs a. 
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Figure 8.- Lateral-directional stability derivatives: nozzles at 60°. 
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and C vs a. 
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Figure 8.- Continued. 
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(a) C £ , C £ , and C £ vs a. 
P r 3 


Figure 9.- Lateral-directional stability derivatives: nozzles at 80° to 90°. 
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(b) C , C , and C vs a. 
n n n 0 


Figure 9.- Continued. 
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(a) Nozzles at 8°. 

Figure 10.- Lateral-control characteristics. 
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(b) Nozzles at 60°. 
Fi gure 10.- Cont inued . 
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A 0.31 TO 0.36 
o 0.38 TO 0.40 
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(c) Nozzles at 80° to 90°. 


Figure 10.- Concluded. 
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Figure 11.— Rolling moment due to the lateral control vs control-wheel angle. 


49 



, per deg 




O 0.1 TO 0.15 

□ 0.21 TO 0.26 

O 0.28 TO 0.31 

PREDICTED 


FILLED SYMBOL - MANEUVER INCLUDED 
CONTROL WHEEL INPUT 



u _4 0 4 8 12 16 20 

a, deg 


(a) Nozzles at 8°. 

Figure 12.- Rudder characteristics. 
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(b) Nozzles at 60°. 
Figure 12.- Continued. 
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(c) Nozzles at 80° to 90°. 
Figure 12.- Concluded. 
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(a) Nozzles angles, 8° and 60°. 

Figure 13.- Comparison of derivatives from the maximum- likelihood and the 

linear-regression solutions. 
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FLAGGED SYMBOL - FIXED VALUE 
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(a) Rolling-moment derivatives. 

Figure 14.- Effect on dispersion of setting fixed values for selected lateral-control 
derivatives: maximum- likelihood solutions; nozzle angles, 80° to 90°. 


55 




C J 

o 0.09 TO 0.14 
O 0.22 TO 0.30 
a 0.31 TO 0.36 
o 0.38 TO 0.40 

FLAGGED SYMBOL - FIXED VALUE 

FILLED SYMBOL - BOTH RUDDER AND CONTROL 
WHEEL INPUTS 

Ci AND C n SELECTED LATERAL 

w w CONTROL 

FREE TO VARY DERIVATIVES FIXED 



a, deg a, deg 


(b) Yawing-moment derivatives. 
Figure 14.- Concluded. 


56 



Ill 



I 


1. Report No. 


2. Government Accession No. 


NASA TP-2033 | 

4. Title and Subtitle 

THE APPLICATION OF PARAMETER ESTIMATION TO FLIGHT MEASUREMENTS 
TO OBTAIN LATERAL-DIRECTIONAL STABILITY DERIVATIVES 
OF AN AUGMENTED JET-FLAP STOL AIRPLANE 


7. Author(s) 


3. Recipient's Catalog No. 


5. Report Date 
January 1983 


6. Performing Organization Code 


8. Performing Organization Report No. 


Jack D. Stephenson 

9. Performing Organization Name and Address 

NASA Ames Research Center 
Moffett Field, Calif. 94035 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


15. Supplementary Notes 


A-8977 

.... 10. Work Unit No. 

T-3478Y 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Paper 

14. Sponsoring Agency Code 
533-02-91 


Point of Contact: Jack D. Stephenson, Ames Research Center, MS 21 1-1, Moffett Field, CA, 94035 

(415)965-6004 or FTS 448-6004 


16. Abstract 


Flight experiments with an augmented jet-flap STOL aircraft provided data from which the 
lateral-directional stability and control derivatives were calculated by applying a linear- 
regression parameter-estimation procedure. The tests, which were conducted with the jet flaps 
set at a 65° deflection, covered a large range of angles of attack and engine power settings. 

The effect of changing the angle of the jet- thrust vector was also investigated. Test results 
are compared with stability derivatives that had been predicted. The roll damping derived from 
the tests was significantly larger than had been predicted, whereas the other derivatives were 
generally in agreement with the predictions. Results obtained using a maximum- likelihood estima- 
tion procedure are compared with those from the linear-regression solutions. 


17. Key Words (Suggested by Author(s)) 
Lateral Stability 
Stability Derivatives 
Augmented Jet Flap Airplane 
Flight Tests 


18. Distribution Statement 

Unclassified - Unlimited 


Subject Category 08 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

Unclassified 

Unclassified 

62 


22. Price* 
A0 4 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1983 


t' 

National Aeronautics and third-class bulk rate 

Space Administration 

Washington, D.C. 

20546 


Official Businef 
Penalty for Pri\ 


V 1 
‘->Cr>T ,)t 
\d 


U,A, 

hid AI,{ 


b 30 

r*0i iCr: 


25 


,n hV:d WdOfMiO.ir 


300903 OS 


SUL ) 


> 



Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 






POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


■ him am i m i in n i i nil l 



